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We study the critical properties of a model of information spreading based on the 
SIS epidemic model. Spreading rates decay with time, as ruled by two parameters, 
e and /, that can be either constant or randomly distributed in the population. The 
spreading dynamics is developed on top of Erdos-Renyi networks. We present the 
mean-field analytical solution of the model in its simplest formulation, and Monte 
Carlo simulations are performed for the more heterogeneous cases. The outcomes 
show that the system undergoes a nonequilibrium phase transition whose critical 
point depends on the parameters e and /. In addition, we conclude that the more 
heterogeneous the population, the more favored the information spreading over the 
network. 
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I. INTRODUCTION 

In the last decades, diverse questions of social dynamics have been tackled by means of 
statistical physics techniques. In fact, simple models allow to simulate and understand real 
problems such as elections, spread of information, vehicle traffic or pedestrian evacuation, 
amongst many others [l| . As a feedback, these issues are attractive to physicists because of 
the occurrence of order-disorder transitions, scaling and universality, among other typical 
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features of physical systems. 

More recently, due to the emergence and popularization of social networks like facebook 
and twitter, as means of information dissemination, there is a growing interest in the study of 
rumor and information spreading in complex networks. For this purpose, a large diversity of 
models have been proposed j3-14]. The majority of these models are based on the standard 
ones of epide mlC spreading like SI, SIS, SIR and their variants Q. The paradtgmattc 
model of rumor spreading is the Daley-Kendal (DK) model [2J , that is conceptually similar 
to the SIR model. The population is divided into three distinct states, namely Spreaders, 
Ignorants and Stiflers. The Spreaders are agents that are transmitting the rumor through the 
population, the Ignorants do not know the rumor and, consequently, they are not spreading, 
and finally the Stiflers are those individuals who know the rumor but have lost interest in 
diffusing it. The transitions between states are given by stochastic rules in the same way as 
in epidemic models. Many extensions of the DK model were studied by the consideration 



of random, scale-free [4-6] and small-world [2J networks for the contact among individuals, 



two different kinds of rumors spreading over the networ 



effects of media 



10], remembering and loss of memory 
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jiiew classes of individuals 



121 ] . impact of human activities 



131 ]. among others. 

Our present motivation is to investigate some of the mechanisms involved in the adoption 
of innovations, new ideas or technologies. This issue may have practica l ap plications such as 
in marketing strategies and have already been target of recent studies [y], Q, We focus 
on the spreading of information of the kind that induces the adoption of a new product 
or idea. Although, for practical purposes, one has primarily in mind the promotion of 
new goods introduced in the market, our proposal may also apply, for instance, to political 
propaganda. For that goal we consider a dynamics of information spreading inspired on the 
SIS epidemic model. Each agent can be in one of two possible states, namely, S (Spreader) or 
R (Restrained). The individuals in state S are those spreading the information through the 
network, whereas the agents in the state R are not transmitting it even if they are aware of 
the information. Notice that for marketing purposes, it is important that the agents become 
enthusiastic transmitters and not only that they know or adopt the technology. Agents 
initially in state R do not have the knowledge. When such an agent is put into contact with 
a new technology, product or ideology, by interaction with Spreaders, the agent becomes 
with rate A an enthusiastic adopter or spreader, trying in turn to convince the neighbors 
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in the network of contacts (friends, relatives, etc.) to adopt the innovation. We denote 
this state as S. Frequently, this enthusiasm is transient and, after some time period, the 
agent spontaneously decays to the state R with rate a, where even though the agent may 
know or use the technology, he/she does not propagate it. Also this state is not permanent 
and the agent could become again a Spreader, however, with a reduced transition rate. To 
take this fact into account, we propose that, once the individual decays to state R, his/her 
"contagion" rate A decreases as A — > e A, with e < 1. It reflects the fact that individuals tend 
to become more and more resistant to spread the information. Nonetheless, the reduction 
of the "infection" rate occurs only a limited number / of times. 

Both parameters, e and I, can be either uniform, i.e., equal for all individuals, or vary from 
one individual to another. In this work we will study the effect of the aging of transition rates 
on the phase diagram of the model. Moreover, we will show the crucial role of heterogeneities 
on the critical behavior and its consequences for information spreading. 

The paper is organized as follows. In Section II we present the general formulation of the 
model and define its microscopic rules. The analytical and numerical results of four distinct 
cases analyzed are discussed in Section III. Section IV contains the conclusions and final 
remarks. 



II. THE MODEL 



We have considered a dynamics of information spreading based on the SIS epidemic 
model, where each agent can be in one of two possible states. In our case, S (spreader) or R 
restrained). The microscopic rules, based on a variant of the SIS model with aging effects 



21] , are the following: 



1. Each individual j in the S state at time t becomes R in the next time step t + 1 with 
probability a; 

2. Each individual j in the R state at time t becomes S in the next time step t + 1 with 
probability Xj(t) if it has at least one neighbor in the S state; 

3. Each individual j starts the dynamics with Xj(t = 0) = Ao- After each transition 
S — >■ R, the spreading probability Xj is updated according to 

X j (t+l) = e j X j (t) , (1) 
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where, for all j, tj < 1 is the factor of reduction of the spreading rate (or probability 
of spreading within a time unit). 

4. In addition, the update given by Eq. ([1]) occurs a maximal number lj of times for each 
individual j. 

Let us recall that, initially, the individuals in the R state do not know the information. 
Actually, this would correspond to a third state (Ignorant agent), but this state only occurs 
at the beginning of the dynamics and it is not attainable later on, as soon as oblivion is 
not taken into account. After a contact with an agent S (spreader), an individual j in state 
R comes to know the information, hence becoming S, with probability Xj (this probability 
is, of course, equal for all individuals at t = 0, i.e., Xj = Xq for all j). After that, the 
individuals will not forget the information, but they can spread it or not. In standard 
models of rumor spreading, the Ignorants do not spread the rumor because they do not 
know it, while the Stiflers know the rumor but they do not transmit it. In comparison with 
the standard states considered in rumor models, the agents in the R state can be identified 
with Ignorants only at the beginning, when they do not have acquired the knowledge yet, 
while they can be identified with Stiflers after each transition S — > R, in which case they do 
have the knowledge. Moreover, concerning the transition rules in standard models, Stiflers 
do not become Spreaders again, unlike in our model. 

Notice, from our rules, that the spreading probabilities are heterogeneous, varying from 
one individual to another. Each agent j who stops spreading the information through the 
network (Restrained) may become a Spreader again, but the probability with which this 
event occurs decreases with time, depending on the agent intrinsic traits, given by (ej,/j). 
After each transition S — >■ R, the spreading probability of a given agent j decreases, up to 
a maximal number of times lj, then, in the next time step, it will be more difficult for the 
contacts of the individual j to "persuade" him/her to spread the information again. 

We have investigated the model on top of an Erdos-Renyi (ER) network with size N 
and Poissonian degree distribution P(k) = e~' fc ' (k) k /k\. In this case, we have computed 
information spreading only in the largest connected component, i.e., the giant cluster. Unless 
otherwise stated, we have considered networks with (k) = 10, for which the probability 



that a given node be 
thermodynamic limit 



ongs to the giant cluster is approximately given by 0.99995 in the 
18]. 
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We will analyze four distinct cases: (i) uniform parameters, i.e., ej = e and lj = I for 
all individuals j; (ii) uniform e and random I; (iii) random e and uniform I, and finally 
(iv) random e and random I. When randomness is considered, the uniform probability 
distribution is used to generate the parameters. 



III. RESULTS 



A. Uniform e and uniform I 



In this case, we have ej = e and lj = I for each individual j. However, there is a certain 
degree of heterogeneity in the system due to the distinct histories of the spreading rates 



Xj. Following the mean- field approach used to treat epidemic models [19], one obtains the 
analytical solution of the model. Initially, let us define S m as the number of Spreaders 
that have performed the transition S — > R exactly m times, where m can take the values 
m — 0,1, Observe that these agents have spreading rates A m = Aoe m . At mean-field 
level, one can assume that after a long time (but before attaining the steady state) all 



individuals will be in either one of two states [20J, namely, either Si or R [21]. Thus, the 



only relevant equation to the time evolution of the system is 

Si = -asi + Xi(k)si (1 - , (2) 

where si = Si/N is the density of Spreaders in the state Si, Xi — Ao e l and we have used the 
normalization condition s + r = 1. Due to the Poisson distribution of ER graphs, in Eq. ([2]), 
we have neglected the fluctuations on the connectivity and made the approximation that 



every node has the same degree 



19] 



In the steady state, si = s and si = s = 0, which leads either to the trivial solution s = 
or to 

S — 1 — 7TTT - 6 ■ (3) 

(fc)A 

This nontrivial solution vanishes at threshold values Ao c in the usual form s ~ Ao — Ao c , 
where the critical points Ao c are given in terms of the parameters by 

V = |y £ -'. (4) 
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FIG. 1: (Color online) Phase diagram of the model, predicted by the mean- field approach, given 
by Eq. in the plane Ao versus e, for (k) = 10 and typical values of I (a). The effect of varying 
(k) is also exhibited for a given value of I, namely I = 5 (b). The region above the curves represents 
the phase where information keeps being spread over the network. 

In other words, for given values of the parameters e and / [22| , there is a critical value of the 
spreading probability Ao c given by Eq. (j3J) separating a phase where the information stops 
being spread (for Ao < AoJ and a phase where a certain fraction of the population remains 
spreading the information (for Ao > AoJ. In Fig. [T] (a) we exhibit the phase diagram of the 
model in the plane Ao versus e for (k) = 10 and typical values of I. One can see that the 
larger I, the smaller the region where the information keeps spreading (the region above the 
curves). This result is easily understood: if we allow the spreading probabilities to decrease 
a large number of times, the probability of the Restrained individuals to become Spreaders 
becomes very small, and it is improbable that the information will remain being spread over 
the network. In fact, this event will occur only if the initial spreading probability Ao is large. 
It is also shown, in Fig. [1] (b), the phase diagram for a fixed value of I (I = 5) and different 
values of the average degree (k). In this case, the spreading phase increases with (k). In 
fact, if each individual has on average a large number of contacts (nearest neighbors), there 
is a greater possibility of spreading the information across the network in comparison with 
the case of a small average degree. 

We confronted the analytical solution of the model with numerical results. We have 
simulated the model on ER random graphs with N = 10 5 nodes and different values of 
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FIG. 2: (Color online) Time evolution of the density of Spreaders, s(t), in the vicinity of the 
transition, for typical values of Ao exhibited in the legend (a). Right at the critical point, the 
density of Spreaders decays as s(t) ~ t , as indicated by the dashed line for Ao = 0.565. These 
results were obtained from simulations of the model with uniform e and I, with (k) = 10, I = 5, 
e = 0.7, N = 10 5 nodes and 100 independent simulations. It is also exhibited the comparative 
phase diagram in the plane Ao versus e for I = 5 and (k) = 10 (b). The squares were obtained from 
the simulations, whereas the full line is the analytical prediction, given by Eq. ([4]). Error bars are 
smaller than symbol size. 



the parameters. For each value of e and I, we have considered different values of Ao and 
100 independent simulations (furnishing configurational averages). As initial condition, 1% 
of the network nodes were set in the S state and the remaining ones in the R state. The 
threshold values A 0c were estimated from the time evolution of the density of Spreaders, 
s(t). Right at the critical point Aq c , this quantity decays in time as the power law s(t) ~ t~ 5 



at mean-field level 



23 



241 ]. As an example, we exhibit in Fig. [2] (a) the time evolution of s 



for 1 = 5, e = 0.7 and some values of Ao in the vicinity of the transition. One can see that 
the above-mentioned power-law behavior can be observed for Ao ~ 0.565. We repeated this 
procedure for other values of e, and we compared the estimated values of the thresholds Ao c 
with the values obtained from the mean- field approach, Eq. We can see from Fig. [2] 
(b) that the considered size (N = 10 5 nodes) gives us a good estimate of the critical points 
Aq c , which confirms that the assumptions made to analytically solve the model are valid. 
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FIG. 3: (Color online) Phase diagram in the plane Ao versus e, for (k) = 10 and / uniformly 
distributed in the integer range [0,10]. The squares were obtained from the simulations, whereas 
the dashed line is just a guide to the eyes. Error bars are smaller than symbols. For comparison, 
the frontier of the uniform case for I = 5 is also exhibited (red full line). 

B. Uniform e and random I 

In this case, ej = e for all j while lj is different for each individual j. In other words, 
we have an additional heterogeneity in the population due to the individual capacity to 
decrease the spreading probability a distinct number of times, i.e., each agent j has a limiting 
parameter lj that is an integer number generated from a uniform distribution in the range 
[0,10]. 

Following the above-discussed procedure, we have analyzed the time evolution of the 
density s(t) of Spreaders for populations of size iV = 10 5 . The critical points Ao c were 
estimated from the power- law behavior s(t) ~ t~ s as in the previous case in Section III. A. 
The critical line in the plane Ao versus e, for I uniformly distributed in [0,10] is exhibited in 
Fig. [3] (points joined by a dashed line). In contrast to what happens in the uniform case 
presented in Section III. A, when I is random the relation between Ao c and e is not a power 
law (see the comparison in Fig. |3]). As a consequence of the heterogeneity of I, one can 
observe that there is a phase transition even for e = 0, which turns the spreading phase 
(above the dashed curve in Fig. [3]) larger than in the case of uniform values of I and e (for 
comparison, we also exhibit in Fig. [3] (full line) the frontier of the uniform case for / = 5, 
corresponding to the mean value of lj). In other words, the transition is not eliminated even 
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FIG. 4: (Color online) Phase diagram in the plane Ao versus /, for (k) = 10 and e uniformly 
distributed in the range [0,1]. The squares were obtained from the simulations, whereas the dashed 
line is just a guide to the eyes. Error bars are smaller than symbols. The inset shows the same 
plot in log-log scale. Fitting data, we obtained the power-law relation Ao c ~ Z 7 , with 7 ~ 0.64 (red 
full line). 

for e = 0. This fact can be understood as follows. In the case of a small value of e, the 
spreading probabilities decrease fast, but there are some individuals for which this decrease 
occurs a small number of times (e.g. for I = 1) or does not occur at all (for 1 = 0). Thus, 
if we consider simulations with a large initial value of A, namely Ao > 0.6, those individuals 
are responsible by the permanent spread of the information through the network. Notice 
that in the limiting case e = 0, the individuals j with lj = (that are around 10% of the 
population) keep their spreading rates equal to the initial value, i.e., they have Xj = Ao at all 
time steps. These individuals can spread permanently the information across the network 
if the initial value of the spreading probability is larger than w 0.6 (see Fig. [3]). In sum, 
for distributions with a given mean value of I, the enhancement of the transmission region 
above the curves is more pronounced if I = is allowed and the larger is the dispersion. 

C. Random e and uniform I 

In this case, we have lj = I while e is different for each individual j. In other words, we 
have an additional heterogeneity in the population due to the individual rate of decrease of 
the spreading probability, i.e., each agent j has a decreasing rate €j that is a real number 
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generated from a uniform distribution in the range [0,1]. 

Again, we have followed the usual procedure and we have analyzed the time evolution of 
the density of spreaders s(t). In the present formulation, as e is random, we have plotted 
the phase diagram in the plane Ao versus I. One can observe in Fig. HJthat the region above 
the curve, for which the information is continuously spread across the network, decreases 
for increasing values of /. In fact, if we increase /, the final spreading rates A; = Ao e l are 
small even if e >> 0, and in this case it is difficult for the Restrained individuals to become 
Spreaders again. If we plot the data in the log-log scale, one case see that the quantities are 
related by the power law Ao c ~ / 7 , with 7 ~ 0.64 (see the inset of Fig. @J. 

D. Random e and random I 

In this case, we have the more heterogeneous instance where both parameters are random: 
e is uniformly distributed in the range [0, 1] and I is an integer number uniformly distributed 
in the range [0, 10]. Although in this case there is no phase diagram to plot, we can discuss 
the criticality of the model at the single transition point Ao c . Our numerical estimate is 
Ao c = 0.31±0.005. In other words, for Ao > 0.31 the information is permanently disseminated 
through the network by a finite fraction of the population. Notice that the threshold is 
relatively small. Thus, the enhanced diversity of spreading probabilities Xj in the population 
favors the propagation of information. 

IV. FINAL REMARKS 

In this work we have studied a model of information spreading on complex networks. 
The model is based on the SIS epidemic model, but the "infection" or spreading rates vary 
with time, decaying with the number of "reinfections". This decay is controlled by two 
parameters, e and /, and we have considered that they can be either uniform or random. 
These features make the population heterogeneous, since the agents may have distinct rates 
of transition between the two possible states, namely, Spreader (S) or Restrained (R). 

We solved the model analytically in its simplest formulation, for constant e and I. In 
this case, the critical spreading rates are given by Ao c = («/ (k)) e" z , where Ao is the initial 
probability (at t = 0) of the transition R — y S, a is the probability of the transition S — > R 
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and (k) is the average degree of the random network. The critical rates define frontiers (for 
different values of I) in the plane A versus e that separate a phase where the information 
stops being spread (for Ao < A J and a phase where a certain fraction of the population 
remains spreading the information (for Ao > Ao c ). This result suggests that networks with 
large average degree favor the spreading of information, as well as populations with a small 
capacity to decrease the spreading rates (i.e., small /), as intuitively expected. All analytical 
results were confirmed by Monte Carlo simulations. 

In the case where one (or both) of the parameters e or I is (are) random, the model was 
analyzed only through numerical simulations. For random / and fixed e, we have observed 
that information spreading is favored and thus the spreading phase, where the information 
is permanently disseminated through the network by a finite fraction of the population, is 
larger than in the uniform case. For random e and fixed I the two phases are of comparable 
size, and the power law A 0c ~ I 1 , with 7 ~ 0.64 arises. Finally, in the case of both parameters 
being heterogeneous in the population, we have found the transition at Ao c — 0.31. Thus, one 
can conclude that the more heterogeneous the population is, more the information spreading 
is favored. In other words, population diversity is an interesting feature to be taken into 
account in models of information/rumor spreading. 

Let us remark that, due to the correspondence between the present model and the SIS 
model with aging, our results can be immediately applied to the latter model. In that case, 
diversity would be malefic, since disease propagation will be favored. 
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